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Abstract 

Studies of charge-charge (ion-ion, ion-electron, and electron-electron) coupling properties for ion 
impurities in an electron gas and for a two component plasma are carried out on the basis of a 
regularized electron-ion potential without short-range Coulomb divergence. This work is moti- 
vated in part by questions arising from recent spectroscopic measurements revealing discrepancies 
with present theoretical descriptions. Many of the current radiative property models for plasmas 
include only single electron-emitter collisions and neglect some or all charge-charge interactions. 
A molecular dynamics simulation of dipole relaxation is proposed here to allow proper account 
of many electron-emitter interactions and all charge-charge couplings. As illustrations, molecular 
dynamics simulations are reported for the cases of a single ion imbedded in an electron plasma 
and for a two-component ion-electron plasma. Ion-ion, electron-ion, and electron-electron coupling 
effects are discussed for hydrogen-like Balmer alpha lines. 

PACS numbers: 52.65.Yy, 52.25.Vy, 05.10.-a 



I. INTRODUCTION 



Theories for the interpretation of spectroscopic measurements in hot plasmas are crucial 
for diagnosing and studying laboratory and astrophysical plasmas. Well known theories such 
as the impact theory Ij or the unified theory Q] for electron broadening will be mentioned 



below. It should be emphasized that in this context a "theory" is a closed ensemble of many 
theoretical developments built on the basis of an initial set of hypotheses, namely a "model" 
(possibly including a few variants). Within each theory, specific analytical and numerical 
methods are assumed for practical calculation of the physical properties that belong to 
that theory. The analysis here addresses the validity of common assumptions regarding 
the many-body Coulomb forces for computing radiative properties. More specifically, the 
lineshapes of ion emitters in fully ionized plasmas are investigated using classical molecular 
dynamics (MD) simulation. The advantage of MD simulation lies in its ability to account 
for all charge-charge interactions without approximation, providing a means to assess their 
relative importance for lineshape calculations. 

Most of the current models for the study of radiative properties of an emitter in a dilute 
medium treat the internal states of the emitter quantum mechanically, but use a semi- 
classical description for their interaction with the environment for the description of the 
environment. The system considered here is the same: a quantum emitter of charge Ze 
imbedded in an equilibrium classical plasma with ions of charge Zie and electrons of charge 
— e at equilibrium. For simplicity, in all of the following it is assumed that the emitter 
couples to the plasma only through its monopole and dipole interactions. 

The line shape is then determined by the Stark broadening of the total electric field at the 
emitter. Three timescales have to be considered together: those characterizing ion motion, 
Tj, electron motion, r e and a typical emitter dipole relaxation time, r^. In all cases Tj is 
much larger than r e due to their large mass ratio. Consequently, the ions do not respond 
directly to individual electron motions but rather to their averaged collective dynamics whose 
primary effect is to screen the ion-ion interactions. Conversely, the free electrons are coupled 
to all other electrons, and a slowly varying ion and emitter charge distribution through 
unscreened potentials. Similarly, the optical electrons bound to the emitter and responsible 
for radiation emission are coupled to the classical charges through a slowly varying screened 
ion field, and rapidly varying electric fields due to the colliding electrons via the dipole 
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interaction. In some current theories the effects of this complex plasma environment are 
represented through the hypothesis of atomic or ionic emitters interacting with ions through 
a Coulomb field screened at the electronic Debye length and by free electrons moving on 
straight trajectories (no electron-emitter or electron-electron interactions) - hereafter called 
independent- electrons. The simplest formulation of the impact theory for electron broadening 
is developed on this basis with three additional constraints for the independent electrons: 1) 
non-overlapping strong collisions (binary collision approximation), 2) positive lower bound 
for the impact parameter range (strong collision cut-off), and 3) a Debye screened Coulomb 
interaction between free and bound emitter electrons. In such theories, collective effects or 
correlations among the various charges are accounted for only in the static Debye screened 
ion-emitter and electron-emitter fields. The objective of our presentation here is to explore 
in a controlled fashion the effects of correlations, both static and dynamic, on properties 
relevant to plasma spectroscopy. 

Lineshape calculations follow two distinct stages: first, the "stochastic" properties of 
the local electric field perturbation at the charged emitter from a classical plasma model 
are determined, and next these properties are used for modeling the dipole emission. The 
first point refers to a determination of charge structure and dynamics, and is essentially a 
nonlinear problem depending on the emitter monopole charge and the plasma parameters - 
a plasma physics problem. The second is connected to the behavior of a quantum emitter 
in a given time dependant stochastic perturbation - an atomic physics problem. The utility 
of MD simulation to treat the plasma physics component "without approximation" is well- 
known j^J, although the application to study correlation effects of electrons with higher Z 
ions is a new direction for hot plasmas 4]. However, in order to assess the importance of such 
effects for spectroscopy, the atomic physics component is essential. The results presented 
here are a first attempt to combine MD simulation and atomic physics for a controlled 
description of both electron and ion broadening j^j]. Earlier applications of MD to spectral 
line broadening were limited to simulating the plasma ion fields j^. Here, we emphasize 
the more difficult task of simulating as well the electrons which couple to the emitter in a 
qualitatively different way due to their attractive interaction. 

The classical plasmas considered for simulation are systems of point charges, either point 
electrons and the emitter monopole in a uniform positive background - hereafter called 
inhomogeneous jellium - or a two component plasma (TCP) of electrons plus the emitter 
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monopoles. It is assumed that quantum mechanical effects are either negligible (as for the 
heavy ions and emitter) or can be incorporated through appropriate modifications of the 
classical Coulomb potential for the electrons. The latter is achieved by using a regularized 
Coulomb potential which remains finite at short distances, accounting simply for the true 
quantum diffraction effects |7[ . This allows derivation of all the desired theoretical properties 
using the laws of classical mechanics, as well as the application of MD. A detailed study 
of the time independent statistical properties of an electron gas with an ion impurity has 
been carried out recently on the basis of such a regularized electron- ion potential Two 
standard classical methods have been used for this study: molecular dynamics (MD) for 
numerical simulation for both dynamical and structural properties, and the hypernetted 
chain approximation (HNC) for the theoretical prediction of structural properties. This 
enables both cross comparisons and interpretations of results. An interpretation of the 
dynamical properties of the electron field autocorrelation function (FCF) at the impurity 
also has been carried out ^| and is considered further here. 

Molecular dynamics is a unique way for providing realistic representations of the stochas- 
tic electric fields determining the dipole interaction of the emitter with its environment. 
Describing these fields by simulation allows one to obtain lineshapes accounting for all the 
correlations between charged particles. Such lineshapes would provide essential reference 
data to benchmark more efficient but phenomenological theoretical models developed for 
plasma diagnostics. Several recent theories for line broadening have been developed for in- 
vestigations at new plasma conditions, including the hot and dense matter conditions needed 
for the presence of highly charged emitters. They suggest that part of the discrepancies found 
with experiments is an inadequate description of the emitter and electron perturber coupling 
mechanisms. The present study provides new information for this discussion. 

The types of questions posed are the following: 1) What is the validity of straight line 
trajectories for the electrons in the presence of a charged emitter (neglect of electron-emitter, 
electron-electron, and electron-ion coupling in the dynamics)? 2) How well are the missing 
correlations of 1) captured by simple screening of the electron emitter interaction? 3) What 
is the importance of multi-electron strong collisions? This study addresses these issues for 
two different classical plasma environments of the emitter. Both are well distinct from those 
resulting from independent electron models assumed for impact theories. The first and 
simplest is the inhomogeneous jellium used for describing electron broadening only. The 
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corresponding study is hereafter referred to as the electron broadening theory (EBT). The 
second is a more realistic two component electron-ion plasma leading to a both ion and 
electron broadening theory (IEBT). The parameter range explored in this work is chosen 
to be compatible with hot and dense plasma diagnostics based on spectroscopy. These are 
conditions similar to those for which there are long standing discrepancies between theory 
and experiment. For a recent discussion and earlier references see reference jsj. 

In section |H] the main points of the classical plasma study for the positive charge emitter 
in the above two plasma environments is described. Section HTTl outlines the dipole relaxation 
stage of the broadening theories reported in this paper. Then, the main results obtained with 
the EBT, i.e. the effects of ion-electron correlations on electron broadening of spectral lines, 
are reported and discussed. Also in Section UTT1 dipole relaxation in a TCP is considered and 
some preliminary results obtained using the the IEBT are shown. The examples reported 
below involve the Balmer alpha line for various hydrogen-like emitters of charge 1, 3, 5 and 
8 and various temperatures. 

II. PLASMA ENVIRONMENT MODELING 

This section provides an illustration of the results that can be obtained from MD sim- 
ulation for the two types of plasma environments considered for the emitter. In both the 
inhomogeneous jellium and the two component electron- ion plasma (TCP), the monopole 
of the emitter is included as a point charge. The dipole operator representing the internal 
states of the emitter is not considered in this section. For the density and temperature 
conditions considered the repulsive interactions of same sign charges effectively excludes 
pair configurations of the order of the DeBroglie wavelength and quantum effects are small. 
Consequently, the ion-ion, ion-emitter, and electron-electron interactions are taken to be 
Coulomb 



where Z\Z 2 is positive. For practical purposes these Coulomb interactions have been 
screened at a distance A = s/2 which is taken to be of the order of the system size s. This 
allows application of the usual periodic boundary conditions in the MD simulation. Since 
interactions among the charges introduce a physical screening at much shorter distances, 
this system size screening does not affect any of the properties being considered here. In 
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contrast, electron-ion and electron-emitter interactions are attractive and the configurations 
within a DeBroglie wavelength become relevant. At such distances the Coulomb interaction 
must be modified within a classical description to account for quantum diffraction effects. 
There are many ways such quantum potentials can be constructed and here we use one of 
the simplest forms [7] 

V ie {r) = -Ze 2 {l-e- r ' s )e- r l x /r (2) 

where 5 = (2irh 2 /m e kBT e Y^ 2 is the thermal De Broglie wavelength. This ion-electron poten- 
tial regularization provides a well defined classical physics for opposite sign charge systems, 
and the various sophisticated classical N-body methods of classical statistical mechanics 
can be applied. The consequences and advantages of such a semi-classical approach for the 
radiative properties of ion emitters will be discussed in the last sections. 

Other parameters of interest are the average electron-electron distance r = (3/47rA^ e ) 1 / 3 
defined in terms of the electron density N e , the electron thermal velocity vq = {2kBT e /m e ) 1 ^ 2 , 
the mean electronic field E = e/r 2 ,, the electron-electron coupling constant T = e 2 /r kBT e , 
and the electron plasma frequency lo p = (ATcN e e 2 /m e ) l l 2 . The coupling between the electrons 
and the impurity is measured by the ion-electron potential at the origin which, in units of 
tq, becomes a = ZT/5. Similar quantities are defined for the ions in the TCP by replacing 
the magnitude of the electron charge and density by the corresponding ion values. 

Molecular dynamics simulations are carried out using a number N of electrons for the 
jellium, or N electrons and N/Zi ions for the TCP. In the first case a point charge Z ^ N 
representing the emitter is included while for the TCP the ions are considered as the emitters 
for simplicity. The simulation is done using periodic boundary conditions on an elementary 
cubic cell of size c large enough to ensure A > \d, where Ad is the Debye wavelength, and 
A < s/2. These conditions assure that the physical screening due to the charge interactions 
occur at a shorter scale than the system size screening, as noted above. Because N must 
be large to agree with the condition Z< JV the impurity- impurity correlation inherent to 
periodic boundary conditions is negligible. Moreover as we are only interested in the charge 
structure and dynamics in the ion vicinity the usual Ewald's sum is not required. Particle 
motion in MD simulation is achieved using a Verlet algorithm with a time step error of order 
(At) 4 : r n+1 = 2r n - r n _x + a„(At) 2 . 
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A. Ion emitter in an electron gas 



This section summarizes the study of a single ion impurity (the monopole of the emitter) 
imbedded in an electron gas. A uniform positive background maintains overall charge neu- 
trality but has no effect on particle motion since, according to periodic boundary conditions 
used in molecular dynamics, the charges move across an infinite system. The motivations 
for doing the impurity case are: 1) it is a well defined and simpler problem than the TCP 
case in comparison, and 2) all the existing theories for electron broadening to compare with 
are based on an impurity approach. The main statistical properties of the electron structure 
close to the ion (e.g., electron density, electric field probability distribution) are derived 
using both MD simulation and the classical hypernetted chain (HNC) integral equations for 
the impurity-electron correlation function. 

Ion-electron and electron-electron coupling affect both the static and dynamic properties 
of the electrons near the impurity ion. The former determines the dominant Z dependence of 
these properties, while the latter controls the screening and other quantitative effects. The 
electron density, electric microfield distribution, and electric field autocorrelation function 
have been studied for various temperatures, densities, and impurity charge number Z. 

Consider first the electron density near the impurity. The simplest prediction is given 
by the Debye-Huckel approximation leading to an enhanced density near the ion and Debye 
screening at large distances. This approximation is expected to fail for increased electron-ion 
coupling (larger values of Z) and/or increased electron-electron coupling. In that case an 
adequate theory is provided by the two component HNC equations [lQj specialized for the 
single impurity case, i.e. for a vanishing ion density. The impurity-electron pair distribution 
function gi e {r) is obtained by numerical solution to these equations as a function of the direct 
electron-electron correlation function c ee (r) calculated separately. Aside from normalization, 
this pair correlation function is the average electron number density, n e (r) = n e g ie (r) at a 
distance r from the ion. Both MD simulation and the HNC approximation are restricted 
to a limited range of values for the impurity charge, density, and temperature. Outside of 
this range frequent electron trapping and non-convergent iterative solution methods for the 
HNC equations occur invalidating the results. Within the parameter space considered here 
MD and HNC results are in good agreement as illustrated in Fig. ^ The enhanced electron 
accretion around the ion impurity is due to the attractive impurity-electron interaction. 
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FIG. 1: Ion-electron pair correlation function at N e = 10 19 cm -3 , T e = 50000i"C, and Z=8. Solid 
line: HNC; circles: MD; dashed: no e-e coupling and Debye screening; dotted: no e-e coupling and 
no screening. 

If this coupling is neglected the electron distribution becomes uniform regardless of the 
electron-electron coupling. Clearly, the independent electron model is a bad approximation 
for the electron density except for Z = 0. Also shown in Fig. Q is the result obtained 
by including the impurity-electron interaction but neglecting electron-electron interactions. 
This misses the screening effects at larger distances due to the electron coupling. For the 
conditions considered, the electron coupling is small, approximately T = 0.116, and the 
simple addition of Debye screening gives quite good results. Other properties involving the 
local electron electric field at the impurity are of primary interest for spectroscopy since 
they determine the dipole coupling to the internal states of the emitter. The electric field 
for the current case is determined from the gradient of the potential in Eq. El For the 
independent electron approximation, the screening of this field is changed to the Debye 
length. Theoretical predictions of the electric microfield distribution can be obtained in 
terms of gi e (r) and the results using HNC have been compared with those from MD, again 
with good agreement. The MD results are illustrated in Fig. El for two values of impurity 
charge Z = 1 and Z = 8. The effect of increased charge number is to decrease and broaden 
the distribution. Also shown is the result for independent electrons. This is the Z = 
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FIG. 2: Field distribution at N e = 10 19 cm~ 3 , T e = 200000^. Solid: Z = 1; dashed: Z = 8; dots: 
independent electrons. 

interaction without electron-electron interactions, but with Debye screened fields. It is seen 
to be similar to the Z = 1 results. As with the average electron density, the dominant 
Z dependence is due to the impurity-electron coupling. However, at fixed Z the electron- 
electron coupling is responsible for important quantitative changes in the location of the 
maximum and width of the distribution. 

The simplest measure of the electric field dynamics is given by the field autocorrelation 
function (FCF). Figure El shows the time dependence for the cases Z = 1 and Z = 8 at 
T e = 200000/^. Also shown is the independent-electron limit obtained for screened fields 
but removing all interactions. As with the field distribution function, the result is similar to 
the Z=l case. The ad hoc screening of the fields assures that the initial value is similar to 
that with electron-electron coupling. This screening also assures that the straight trajectory 
dynamics is effectively restricted to a Debye sphere. However, the electron-impurity coupling 
changes this agreement dramatically at larger values of Z. Qualitatively, the behavior of the 
correlation function can be summarized as follows: For increasing impurity charge number 
the initial value < E 2 > /Eq increases while the correlation time decreases. Both effects 
are missed by the independent electron approximation, but are regained qualitatively if only 
impurity-electron coupling is retained (without electron-electron interactions). However 
electron-electron coupling is required for quantitative effects as indicated below. 
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FIG. 3: Field autocorrelation function, same as Fig. 

Now consider the effects of electron-electron coupling in more detail for the case Z = 1. 
Figure E] shows the results for two different temperatures, T e = 2000007^ and T e = 50000fT. 
This corresponds to two different values for the electron-electron coupling constant, V = 
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FIG. 4: FCF for two temperatures at Z=l. Dash: no coupling; solid: all couplings. 

0.029 and T = 0.116, and two different values for the electron-ion coupling constant, a = 0.5 
and a = 1, for the high and low temperature cases, respectively. Also shown for comparison 
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are the results for the independent electron model. The initial value of the FCF is seen 
to increase at constant Z for increasing temperature. This can be understood since the 
regularized potential leads to a covariance proportional to 1/5 ~ T 1 / 2 . The correlation time 
decreases at higher temperatures, in part at least due to the increased thermal velocity. 
As expected, the deviation of the independent electron model from these results increases 
with increased electron coupling. The apparent agreement for the high temperature case 
is misleading in this presentation, since the percent discrepancies are large over the entire 
time interval as will be shown in Fig. [7| below. The effect of electron-electron coupling 
on the dynamics is therefore important. This is more evident for the lower temperature 
comparison. 

The strong electron-impurity coupling case of Z = 8 is shown in more detail in Fig. A 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

time (ra 1 ) 

FIG. 5: FCF anti-correlation for Z = 8 and T e = 10000CLFC 

striking new feature is the anticorrelation (negative FCF) at intermediate times. This effect 
increases for larger Z and also for decreasing temperature, corresponding to the enhanced 
electron-ion interactions. This increases the probability of electron configurations near the 
ion. Many of these lead to temporary orbiting trajectories with associated large oscillatory 
electric fields. When such trajectories dominate the contribution to the FCF, their changing 
fields on average will lead to an anticorrelation as observed in Fig. |S| The occurrence of such 
events is easily seen in the MD simulations, as shown in Fig. |Hl The total field is dominated 
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by the oscillatory fields of a few close orbits. These oscillations persist on the time scale of the 
FCF, but are seen to break up on a time scale an order of magnitude longer, due to electron- 
electron interactions. Thus, the effect of anticorrelation is due to quasi-bound orbits, but 
their correct equilibrium statistical sampling requires the electron-electron interactions to 
assure the properly weighted ensemble of temporary orbiting trajectories. 
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FIG. 6: Example of temporarily oscillating electric field at the ion. Z = 8 and T e = lOOOOOi^. 



A general remark can be addressed about the MD approach in the context of this dis- 
cussion. The field history sampling is started after a thermal equilibration stage intended 
to reach a stationary energy state. Independent samplings are obtained starting from a new 
uniform random configuration of particles followed by a thermal equilibration stage. This 
process is repeated many times in order to build a very large sample set. If the electron- 
electron coupling is neglected, the equilibration stage does not favor the electron accretion 
around the ion impurity in the correct way. As a result the sample set obtained is certainly 
not the same as the sample set obtained when the electron-electron interaction is accounted 
for. An improved sampling could be built using initial configurations chosen in order to 
agree with the nonlinear electron distribution around a positive charge. This way has been 
followed to develop a simple theoretical model for the FCF ^]. 

It will be seen later that the spectral line width is proportional to the time integral of 
the FCF. It is difficult to predict the Z dependence of this integral a priori due to the above 
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competing effects of increasing initial value and decreasing correlation time. Consider again 
the FCF for Z — 1 in Fig. 0] at two different temperatures. The time integral of these 
correlation functions and another at lOOOOOfC are shown as a function of temperature in 
Fig. [7| Also shown are the corresponding results for the independent particle model. It is 
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FIG. 7: FCF integral versus temperature for Z = 1. white triangles: all couplings, black triangles: 
independent-electrons. 

seen that the FCF integral decreases with temperature, so evidently the dynamical effect 
of a decreasing correlation time dominates the increasing initial correlations. The effects 
of electron-impurity and electron-electron coupling leads to an approximately 15 percent 
decrease in the integral at all temperatures relative to that for the independent electron 
model. 

Finally, the Z dependence of the FCF integral is demonstrated in Fig. |HJ Again dynamical 
effects of a decreasing correlation time and anticorrelation dominate, leading to an approxi- 
mately linear decrease of the integral with increasing Z. Consider first the value Z = (no 
impurity-electron coupling). The triangle represents the results with electron-electron cou- 
pling, while the line is the independent electron model. It is seen that the use of the screened 
field in the latter accounts well for the relevant electron-electron coupling. This is consistent 
with the discussion following Fig. El Also shown is the independent electron model without 
using the screened fields (open circle). This shows clearly that electron-electron coupling is 
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FIG. 8: FCF integral versus emitter charge at T e = 200000^. Squares: all couplings; dashed: 
independent-electrons; triangle: all couplings with Z = 0; circle: unscreened uncoupled electrons 
with Z = 0. 

almost a 10 percent effect at Z = 0. This is significant since the electron-electron coupling 
constant is quite small (r = 0.029). For Z — 1 the electron-emitter interaction is quite im- 
portant, decreasing the integral by almost 15 percent below the independent electron value. 
This coupling continues to dominate for increasing Z. As noted above, for Z > there is 
an additional importance of the electron-electron coupling in the MD which is needed to 
establish the correct weight for the increasing number of quasi-bound orbits. 

B. Ion-electron two component plasma MD simulations 

In this subsection MD studies of an equilibrium, charge neutral two-component plasmas 
(TCP) of iV electrons and N/Zi ions are reported. Charge neutrality is guaranteed by 
the charge concentration. Preliminary classical TCP MD simulations in this context of 
applications to spectroscopy have been reported elsewhere The challenging aspect of 
such simulations is to move slow and fast particles at the same time. Such a study of fast 
and slow processes requires simulations stable over a long period of time, long enough to 
account for ion motion and based on a time step compatible with electron motion. This 
problem has been addressed mostly in the context of hydrogenic plasmas [11] and is extended 
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here to higher Z ions. The emitter is taken to be one of the ions of the TCP. 

It is worth noting at this point the distinction between classical MD applied to the TCP, 
where both ions and electrons move dynamically, and "quantum MD" where the ions move 



classically on a self-consistent potential energy surface of the electrons 3, 13 1- In the latter, 
at each time step for the ions the ground state configuration for the electrons is determined by 
some quantum method (e.g., density functional theory). Subsequently the forces between 
the ions in this electronic configuration are calculated and the ions are moved along the 
corresponding classical trajectory. The process is repeated at the next time step. The two 
classes of simulations appear complementary, with quantum MD being more appropriate for 
lower temperatures and semi-classical MD more appropriate at higher temperatures. 

We first verify that the MD simulation properly describes the expected ion-electron 
screening mechanism. For this purpose the simulation of a neutral mixture of hydrogen- 
like carbon ions Z = 5 in electrons has been carried out with the temperature T e = 50000fC 
and the electron density N e = 10 19 cm~ 3 . As discussed following Eq. a large cutoff length 
A comparable to the system size has been used in the potential to allow the periodic bound- 
ary conditions. Since the physical screening resulting from the collective effects of these 
"Coulomb" interactions is at a much smaller scale, this cutoff does not affect the study of 
that screening. A comparison of the ion-ion pair distribution function for the TCP with 
"Coulomb" interactions with that for a corresponding ion OCP using Debye screened inter- 
action is made. In the OCP system there are no electrons and ion-ion Coulomb potential is 
screened at the Debye length phenomenologically, while such screening must be generated 
directly by the electrons in the TCP simulations. Figure |H1 shows the rather good agreement 
between the two simulations, giving an indirect evidence of the effectiveness of the ion-ion 
dynamic screening by the electrons (also shown is the OCP without screening to quantify 
the effect being studied). The TCP result shows a small enhancement of the screening 
mechanism over that of the OCP at shorter distances. This behavior suggests that the short 
range ion-ion interaction could be well described using an effective ion charge Z < Z. An 
rough estimate not shown gives Z ~ 4.7. 

Evidence of screening also can be demonstrated via the dynamics of the total field due 
to both ions and electrons at one of the ions (considered to be the emitter). MD simulation 
allows determination of statistical data for each particle type involved. In the electron-ion 
TCP two kinds of time dependencies coexist, the high frequency and the low frequency 
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FIG. 9: Pair correlation function. Solid: OCP; circles: TCP; dots: no screening. 

dynamics related to electron motion and ion motion, respectively. Both of them are of 
interest for probing or discussing the simulation procedure accuracy. For the high frequency 
data the electron field autocorrelation function at the ion has been calculated. This is 
similar to the property calculated in the subsection above for jellium, except that now it is 
in the presence of many other ions. Figure El shows a comparison of the field autocorrelation 
function for the field seen by a single ion in electrons (jellium) and for the TCP. Despite some 
noise due to insufficient sampling it can be observed that for the TCP the field correlation 
function does not go to as for the impurity in jellium case. The timescale used in this 
simulation is the inverse electron plasma frequency, compatible with the high frequency 
mechanisms of the electrons. On this time scale the electron FCF tends to a constant value 
rather than zero, denoting the occurrence of a low-frequency component in the electron field, 
due to its coupling to the slower ion dynamics. Contrary to the impurity ion in jellium, a 
given ion is no longer a center of symmetry for the average electron structure. Rather, 
this average electron structure is a blurred image of the random ion configuration which 
remains quasi-static at this time scale. This low frequency component disappears, as shown 
on Fig. ^2 if the ions are located on the nodes of a cubic array where each ion is a center 
of symmetry of the infinite system. As expected, on the time scale of this low-frequency 
dynamics this seemingly constant value goes to 0. 

In many theoretical formulations of lineshapes, such as the fast fluctuation approximation 
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FIG. 10: Electron field autocorrelation function. Dashed: impurity; circles: TCP; solid: cubic 
array. 

of the next section, it is assumed that the electron electric field relaxes quickly relative to 
the time scale of the emission or absorption process. One consequence of the occurrence of 
a low-frequency component in the TCP electron field is that these approximations are no 
longer valid. 

The FCFs generated by a one component plasma of ions interacting through a Coulomb 
potential 1) with a screening length Xd = Debye length, 2) with a box screening length A > 
A d ar e compared with the FCF of the low frequency field component of the ion plus electron 
total field on Fig. EI The good agreement of the two curves and the large discrepancy with 
the no screening result can be interpreted again as the proof that the screening mechanisms 
take place properly in the TCP molecular dynamics simulations. 



III. ELECTRON-ION COUPLING EFFECTS ON LINESHAPE 



A. Dipole relaxation mechanisms 

Linear response methods can be applied to the dipole relaxation of a quantum emitter in 
a classical plasma environment Q|. According to the energy involved it is postulated that 
the perturbation by the environment never induces transitions between the two manifolds 
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FIG. 11: Total field autocorrelation function. Dots: TCP simulation; solid: ion OCP simulation; 
dashed: ion OCP no screening. 

for which the transition is considered within the characteristic time of the emission - no 
quenching approximation. This approximation therefore neglects emitter-bath coupling via 
the radiated dipole transition and any contribution of the dipole transition moment at 
equilibrium. As a result, the spectral profiles of dipole spontaneous-emission are given in 
terms of the autocorrelation function C(t) of the system dipole moment q(t), via a Fourier 
transformation: 

1 f°° 

I(u) = Re- / dt C{t) e luJt , (3) 

7T Jo 

C(t) = ti {ir(t)pqU(t).q} (4) 

The trace is taken over an equilibrium Gibbs statistical ensemble p for the total system of 
plasma plus emitter. The Hamiltonian is taken to be the sum of the plasma Hamiltonian 
(including the point monopole of the emitter) H p , the Hamiltonian for the internal states of 
the emitter H e , and their interaction due to the plasma electric field coupling to the internal 
emitter dipole. For simplicity here, the Doppler broadening due to emitter center of mass 
motion has been assumed to be statistically independent of this dipole broadening. The 
operator q(t) = U(t)qU*(t) obeys the Heisemberg equation 

= [H e - E(f) • d, q(*)] (5) 
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The second term in the commutator describes the perturbation due to the total electric field 
of the plasma E(i). Due to the no quenching approximation, d is restricted here to the part 
of the internal emitter dipole operator coupled to this field. The time dependence of E(£) 
is generated by the plasma Hamiltonian alone. This shows most directly the utility of MD 
simulation in addressing the most difficult part of a line broadening problem: identifying 
properly and completely the environment of the emitter. In the present case this is given 
by an electric field history. To represent an equilibrium ensemble for the environment a 
collection of field histories for different initial conditions is given, the Schrodinger equation 
solved in each case, and a simple algebraic average is performed over all such solutions to 
get the line profile. In practice the numerical solution to the Schrodinger equation can be 
straightforward or complicated depending on the electric field history. This is the case for 
high Z radiators for which the quasi bound orbits described above can give rise to high 
frequency, strong fields (e.g., as in Fig. EJ. The calculations require a fast integration 
process based on a polynomial development using the so called Euler-Rodrigues coefficients. 
A detailed presentation of the method applied to hydrogenic lines has been given elsewhere 

0. 

The existing theoretical analysis generally entails further approximations not necessary 
in the above MD/Schrodinger equation simulation. Due to the small electron-ion mass ratio, 
Stark effects resulting from slow ion and fast electron micro-fields are generally considered 
separately. The Stark effect due to ions is approximated by a static electric field, while 
electron broadening is described as a dynamical process resulting from the high frequency 
fluctuation of the electron electric field. The usual approach for electron broadening [14| for 
ion emitters in plasmas is a binary collision model involving independent electrons moving 
on straight trajectories at constant velocity. Electron correlations are accounted for only 
indirectly by screening the electron-emitter potential at the Debye length. There are other 
theoretical models with different approximations, but for the discussion here we will primar- 
ily compare simulations for two cases: 1) the electrons are considered as free particles that 
do not interact with each other but perturb the ion emitters with a field obtained from the 
potential written in Eq. (J2J) with A = A^; and 2) the electrons have Coulomb interactions 
with each other and couple to the emitter with a field obtained using Eq. (J2J) with A = s/2. 
The latter, of course, is the correct treatment of correlations. 

A more compact formulation of the problem is obtained using a Liouville operator repre- 
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sentation Q| allowing to disentangle the evolution operator from the internal emitter dipole 
operator q. Liouville operators are defined according to the following relation: Lq = hi^o, •]■ 
In this representation the dipole correlation function can be rewritten C(t) = Tr(q-U(t)pq). 
Equation El is replaced by the following stochastic equation for the evolution operator alone: 

= (L + D • E a (t))U a (t), U a (0) = 1 
1 - 

U(t) = {u a (t)} av = -J2 u f«® ( 6 ) 

0=1 

where E a (i) is one of the electric field histories, D is the Liouville representation of d and 
U(t) denotes the solution of the stochastic equation. 

A useful approximation described in Appendix A, the fast fluctuation limit, can be ob- 
tained when the characteristic time for the field correlation function is much smaller than 
the typical relaxation time of the physical process investigated - here the dipole -17]. In this 
limit the stochastic equation reduces to: 

^1 = (L - W)U(t) (7) 

where W is proportional to the field correlation function integral / °°{E Q ,(0)-E Q ,(t)} Q <it. Thus 
questions regarding the effects of charge-charge coupling on the linewidth can be addressed 
directly from the simulation of the field autocorrelation function as described in the previous 
sections. 



B. Electron broadening 

The purpose of this section is to apply the EBT developed herein for discussing the 
occurrence of charge-charge coupling effects in electron broadening. This issue is motivated 
in part by recent experimental and theoretical studies performed on a series of lines emitted 
by ions of increasing charge P|. The specificity of these lines relies on the absence of ion 
broadening. However, in contrast, the present study concerns lines with ion broadening so 
the generalization of the conclusions here to all kinds of lines would be inappropriate. Ion 
lines with non negligible ion broadening have been selected as they fit also the objective of 
studying Stark broadening in two component plasmas accounting for complete charge-charge 
coupling. 
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According to the EBT, the study of electron broadening relies first on the construction of 
large field sample-sets at the emitter generated by molecular dynamics. Then, the lineshape 
or simply the full width at half maximum (FWHM), j(Z,T e ), is obtained after numerical 
solutions of the Schrodinger equation and its average. Due to the various characteristic 
times of the problem, a given field history can be split into a number of independent samples 
depending on the final purpose. In the present case, the dipole relaxation simulation requires 
samples much more longer than those necessary to obtain the field statistical properties. In 
addition, the study of electron broadening as a function of the emitter charge can be limited 
by calculation costs, since the dipole response of the emitter decreases as the inverse of the 
square of the emitter charge. For instance, to carry out an equivalent study with respect 
to the noise inherent to this kind of statistical average for emitters of net charge Z — 1 
and Z — 8, it would be necessary to generate a field sample set and then integrate the 
Schrodinger equation over a time about 20 times longer for Z = 8 than for case Z = 1. 
Accordingly, the present study is limited to cases relevant for the fast fluctuation limit 
(FFL) allowing a simpler theoretical representation of the line width and thereby avoiding 
the full integration process. It should be noted that, for a given temperature, if conditions 
for the FFL are fulfilled for a chosen line in the case Z = 1 they are even more relevant 
for the same line in the case Z = 8. However, the dipole autocorrelation will be calculated 
by complete integration of the Schrodinger equation in a few cases in order to check the 
FFL and for cross comparisons. In order to probe the statistical properties of the local 
microfield, Balmer-alpha lines for hydrogen like emitters of charge Z = 1, 3, 5, 8 have been 
chosen together with the three plasma temperatures T = 5 x 10 4 , 10 5 , 2 x 10 5 K selected 
for this study. Both the dipole moments of the lower and the upper level manifolds are 
coupled to the perturbing field giving rise to the relaxation process and line broadening of 
the emitted dipole transition n = 3 — * n = 2. Qualitatively, for the Balmer alpha line, 
the electric field perturbation induces two opposite broadening mechanisms, i.e. the line 
broadening increases with the covariance characteristic of the field strength while, due to 
some time averaging inherent to the dipole relaxation process, it decreases with an increasing 
fluctuation dynamics of the perturbation. 

Comparisons of the FWHM are reported in Tab. HI where ^ffl is the fast fluctuation 
limit result and jse results from simulation i.e. from a full integration of the Schrodinger 
equation. The discrepancy is smaller than 10% of the FWHM. In addition the exponential 
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TABLE I: Balmer alpha FWHM (s 1 ) versus temperature for Z = 1. 

Temperature Jffl Ise 

5 x 10 4 K 10.7 x 10 12 11.8 x 10 12 

10 5 K 9.88 x 10 12 10.2 x 10 12 

2 x 10 5 K 9.32 x 10 12 9.93 x 10 12 




FIG. 12: Balmer-alpha FWHM, Z = 1, T e = 200000 K. Squares: independent-electrons; solid 
line: independent-electron linear fit; circles: charge-charge coupling; black triangle: complete 
integration; dots: 1/(Z + 1) 2 line. 

decay of C(t) predicted by the FFL also is confirmed for the simulations (not shown). 

For the high temperature case, these results are plotted on a log scale in Fig. EI versus 
the nucleus charge, together with the corresponding independent electron calculations (free 
electrons moving on straight trajectories that perturb the emitter via the potential defined 
in Eq. |21 screened at the Debye length). The field sampling generation cost by the MD 
technique is quite low for noninteracting particles. Therefore, these calculations can be 
performed by integration of the stochastic equation with a satisfactory noise level. The 
FWHM calculations fit the 1/(Z + l) 2 behavior even without coupling to the emitter due to 
the atomic physics dependence (dipole matrix elements) of the linewidth. These calculations 
confirm the capability of the simulation technique to provide relevant results even in the 
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fast fluctuation limit domain. The entire discrepancy between the two results of Fig. ^] 
is therefore a consequence of the charge-charge coupling in the time integral of the FCF 
as a function of the emitter charge (i.e., the effect discussed in Fig. [SJ. This discrepancy 
increases with decreasing temperature as shown in Fig. [7| 

C. Coupled ion-electron broadening 

It has been shown in a previous section that classical two component MD simulations 
provide relevant results regarding charge structure and dynamics. The plasma conditions 
have been selected to emphasize screening mechanisms in order to demonstrate significative 
effects. Other results obtained for weaker coupling conditions lead to the same observations. 
These investigations ensure that ion plus electron field samplings are properly generated by 
the MD simulation for the plasma conditions of this work. Like for the EBT, field sampling 
is the first stage of the ion electron broadening theory (IEBT). In the second stage these 
fields are used for integrating the Schrodinger equation in order to synthesize lineshapes. 
This protocol requires a numerical cost much smaller than for the impurity case because 1) 
a number of samples are generated at the same time using each ion as the emitter, and 2) 
the integration time required for C(t) is shorter when ion broadening is accounted for. Two 
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FIG. 13: Helium hydrogen- like Balmer alpha line T e = 50000 K. 
field data sets are used in order to look for evidences of ion-electron coupling mechanisms at 
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the lineshape level. The two data-sets are built using TCP molecular dynamics techniques, 
1) with all the charge-charge interactions accounted for and, 2) taking into account only the 
ion-ion interactions assuming a Debye screening for the ion-ion potential and free electrons 
moving on straight trajectories. The observed discrepancy in Fig. EI] is a decrease of 6 — 7% 
in the full width at half maximum when both the ion-electron and the electron-electron 
coupling are accounted for. 

IV. CONCLUSION 

Semi-classical models involving a quantum emitter perturbed by a stochastic classical 
electric field have been developed to calculate emitter dipole relaxation. On this basis 
charge correlation effects relevant for plasma spectroscopy have been explored. Such ref- 
erence results from ideal numerical experiments are indispensable, particularly when the 
plasma conditions are such that the usual theoretical approximations appear questionable. 
The calculations reported here are for two distinct models, with and without ion-electron 
and electron-electron couplings. They give rise to comparisons of two sets of line broaden- 
ing calculations. The observed residual discrepancies confirm that effects of charge-charge 
couplings on lineshapes can be non negligible. 

Common electron broadening theories do not account for ion-electron couplings. The 
main reason is that for most of the ionic lines (normal lines) electron broadening is a small 
effect with respect to ion broadening. This has not motivated the development of theories 
accounting for the complex problem of the quantum properties of a charged emitter under- 
going a perturbation resulting from the electrons around the charge. However, this interest 
exists for lines whose broadening is dominated by the electrons. In this work charge-charge 
coupling effects on ionic lines have been investigated considering electron broadening of lines 
with ion broadening. It is worth analyzing the main results obtained here with the impact 
theory for comparison. This theory is founded on the hypothesis that strong collisions do 
not overlap in time. Collision operators for electron broadening have been derived on this 
basis using either linear or parabolic electron trajectories [18] . For the electrons attracted to 
the ion impurity the impact theory postulate of binary encounters is no longer correct as the 
close electron trajectories depend on the other close electrons via kinetic energy exchanges, 
particularly in the case of temporary orbiting trajectories. Accordingly it can be concluded 
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that ion-electron couplings cannot be treated ignoring electron-electron couplings. In terms 
of the electric field dynamics, these couplings induce an increase of both the covariance and 
the field fluctuation rate when the emitter charge or the temperature increase. These effects 
result in a decrease of the width of the studied lines in comparison with the uncoupled case. 

Both the EBT and the IEBT have been designed to account for all charge-charge cou- 
plings. They are based on the numerical integration of the Schrodinger equation over a large 
electric field sample set generated by molecular dynamics techniques. For the selected lines 
and plasma conditions the fast fluctuation limit provides a simple theoretical expression for 
electron broadening in terms of the time integral of the field autocorrelation function. This 
approximation has been checked by comparisons with a few results obtained by full integra- 
tion of the Schrodinger equation. In parallel to this work MD simulations and theoretical 
studies of the FCF for an ion impurity in jellium have been carried out with preliminary 
results reported in |8] . Thus, the FCF and the fast fluctuation limit provide an efficient pre- 
diction of the linewidth for this simple impurity case. In contrast, it appears that integration 
of the Schrodinger equation appears easier when both ion and electron broadening occur to- 
gether, using two component plasma MD simulations for the field histories. It has been 
shown here that such TCP simulations reliably yield the proper dynamic electron screening 
of the ion-ion forces. 

Finally, it can be noted that the direct integration of the Schrodinger equation for line- 
shapes can be rather intricate and expensive. Despite the apparent difficulty of this step by 
step integration process, it appears robust enough to give accurate results that can be seen 
as ideal experiments providing reference data appropriate to discuss or benchmark other 
theories. 

Acknowledgments 

Support for this research has been provided in part by an European action integrated 
Picasso plan. The Spanish group has been financed by de DGICYT (project FTN2001- 
1827) and the Junta de Castilla y Leon (project VA009/03). The research of J. Dufty was 
supported by U. S. Department of Energy Grant DE-FG02ER54677 and by the University 
of Provence. 



25 



APPENDIX A: FAST FLUCTUATION LIMIT 



In this Appendix an outline of the fast fluctuation limit for the line profile is given. The 
terminology "fast fluctuation limit" used instead of the alternative "impact limit" has been 
preferred to avoid confusing it with the impact theories. The later relies on the calculation of 
an average binary collision - the collision operator. The impact theories generally postulate 
non overlapping strong collisions, which is not a limitation of the fast fluctuation limit 
obtained here. 



1. Characteristic time scales 



Three characteristic time scales have to be considered. First, the typical correlation time 
of the perturbing field E(t), and in consequence, the characteristic time of L(t). This time 
scale is ruled by the kinetics of the charged particles in the plasma. Its order of magnitude 
is 

re ~ - (Al) 
^0 

where r denotes typical inter-particle distance, v mean thermal velocity. 

A second time scale is fixed by the correlations of the dipole-moment q(£). The spectral 
width is determined by this lifetime, r^, which is, in a way, the "unknown variable" of the 
problem. 

Finally, a time scale is fixed by characteristic values of L : 

— ~L~ e £ £ (A2) 

do is the Bohr radius, n the principal quantum number of the upper group of states. Evo- 
lution of the dipole moment q(£) would be fixed by this frequency scale if the perturbing 
electric field were stationary. 

The relationship between th and r c will determine the relevant physical phenomenon in 
spectral line broadening. r c > t h is the condition for "quasistatic" broadening which is a 
close limit for ion broadening. Both shape and width of the spectral line are fixed by the 
statistical distribution of the perturbing fields, then ~ Ta- 
li r c <C t h the evolution of the dipole moment q(t) is much slower than the evolution of 
the perturbing fields. This condition fits the electron broadening mechanisms. In this case, 
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perturbations are less efficient, since the emitter responds to a time average of the electric 
field smaller than its statistical typical value. It is the fast fluctuation regime. In the next 
section it will be seen that in this regime the dipole autocorrelation fits well a decreasing 
exponential whose lifetime ^> t# > r c , is determined by the autocorrelation function 
of the perturbing field. It should be noted that in the present study the ratio r c /r^ varies 
from a few hundreds to a few thousands. Such a large ratio guarantees that the FFL is 
totally relevant as the discrepancy of the dipole autocorrelation function to an exponential 
(or equivalently the discrepancy of the lineshape to a Lorentzian) is non significant. 

2. Fast fluctuation limit 

The dipole correlation function of the text can be rewritten as 

C(t) = tr { q(i) • qp } = tr e ( tr p { e Lt pq }) • q (A3) 

The trace has been separated into that for the emitter and plasma subspaces, stationarity 
of the Gibbs ensemble under the dynamics has been used, and the Liouville operator L 
generating the dynamics has been introduced. It is the sum of the Liouville operator for 
the emitter in the relative coordinate system, L e , the Liouville operator for the plasma 
including the center of mass (point monopole) of the emitter, L p , and the coupling between 
the internal emitter states and the plasma, Li 

L = L e + L p + L t (A4) 

Define the interaction representation generator by 

e Lt = e {L e +L p )t u( ^ = e (L e+ L p )t rp exp f* dt > L .^ ( A5 ) 

Jo 

Here T is the time ordering operator with largest times to the left. Then 

C(t) = tr e ( e^ +L ^ e 6 « / e q ) • q ( A6) 

where 0(t) is the average generator for the interaction dynamics in the atomic subspace 

e 6 « = tr p T exp f dt' L^pf; 1 = (T exp f dt'L^t'))^ f e = tv p p (A7) 
Jo Jo 

The leading terms in a cumulant expansion of this average are 

0(t) = 6 (1) (t) + (2) (t) + .. (A8) 
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with 

0«(t)= f dt\U(t')) p (A9) 
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d^(t) = f dt' f dt" (L,,(t')W)> P = t fdt'fl-^) {L % {t)U{t')) p (A10) 
Jo Jo Jo V 
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where 

Zi(0 = Li(0 - tr p L^)pf: x (All) 

The fast fluctuation limit corresponds to the case where the perturber dynamics varies 
rapidly on the time scale of the dipole correlation function C(t). Since the average is per- 
formed only over the plasma degrees of freedom this time scale is controlled by the perturbers. 
Then the fluctuation in Eq. (|A10|) decays to zero for if = r <C i and the factor t'/t in the 
integrand can be neglected 

d^{t)^t [ t dt'(L i (t)L i (t , )) p (A12) 
Jo 

For similar reasons, higher order terms in the cumulant expansion are expected to be of 
higher order in r/t and therefore negligible. 

Further simplifications occur with additional assumptions. If the interaction is only 
through a coupling of the emitter dipole to the plasma field then L, = — E • D. Also, if the 
emitter-plasma interactions are neglected in pf~ l then 

Pie'' - Pp (A13) 

and the fluctuations become 

(Li(t)) p - 0, 

(Zi(Oli(O)p - ~r(t-i')D(i)-D(i') (A14) 

where T(t) is the electric field autocorrelation function for the plasma alone 

r(f) = (E(t) ■ E(0)) p (A15) 

and the dipole operator in D(t) has a time dependence due to the emitter alone. Finally, in 
the case of an emitter with upper and lower degenerate manifolds such as the Balmer alpha 
line without fine structure this time dependence of the dipole operators can be neglected as 
well. The result quoted in the text is then obtained 

C{t) = tr a (e^VU" dt ' V{t,) D D ) / e q) • q ( A16) 
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